Singular features in noise-induced transport with dry friction 
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We present an exactly solvable nonlinear model for the directed motion of an object due to zero- 
mean fluctuations on a uniform featureless surface. Directed motion results from the effect of dry 
(Coulombic) friction coupled to asymmetric surface vibrations with Poissonian shot noise statistics. 
We find that the transport of the object exhibits striking non-monotonic and singular features: 
transport actually improves for increasing dry friction up to a critical dry friction strength A* and 
undergoes a transition to a unidirectional mode of motion at A*. This transition is indicated by a 
cusp singularity in the mean velocity of the object. Moreover, the stationary velocity distribution 
also contains singular features, such as a discontinuity and a delta peak at zero velocity. Our results 
highlight that dissipation can in fact enhance transport, which might be exploited in artificial small 
scale systems. 
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I. INTRODUCTION 

Many transport processes in physics and biology rely 
on the generation of directed motion in the absence of 
an externally imposed force gradient. In recent years 
much research effort has been devoted to understanding 
the physical principles underlying this kind of fluctua- 
tion induced motion, and to the question of how one can 
implement these principles in artificial systems on micro- 
and nanoscales PHI] . Unidirectional motion results quite 
generically from the combination of (i) non- equilibrium 
energy input, e.g., due to chemical or mechanical fluc- 
tuations, and (ii) a spatial or dynamical asymmetry. In 
the paradigmatic example of a Brownian ratchet [6] , 
thermal energy alone is not sufficient to generate a par- 
ticle drift in a ratchet shaped potential due to the re- 
striction imposed by the second law of thermodynam- 
ics, as was discussed already in the classical works by 
Smoluchowski, Feynman, and Huxley (see, e.g., [3] and 
references therein). A non-zero drift results only when 
detailed balance is broken, which can be induced, e.g., 
by periodic switching between low and high temperature 
states, such that the Brownian particles can diffuse more 
easily over the potential barriers causing a net motion in 
a direction prescribed by the asymmetry. 

In this letter we consider a transport model in which 
directed motion of an object is generated by coupling 
noisy asymmetric vibrations, taken specifically as Pois- 
sonian shot noise (PSN), with dry (or Coulombic) fric- 
tion. Both asymmetry and non-equilibrium energy input 
are in this case due to the vibrations, but by itself this 
would not be sufficient to generate directed motion. In 
addition, the moving object has to exhibit a nonlinear 
resistance to the imposed force, leading to a hysteresis in 
the contact line that rectifies the asymmetric vibrations. 
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This nonlinear response arises due to the dry friction be- 
tween the object and the surface. 

Our model thus combines two rather generic ingredi- 
ents: on the one hand, PSN represents stochastic im- 
pulses that occur with a certain frequency and can be 
considered as a generalization of the usual Gaussian 
noise, to which it converges in an appropriate limit [7]. 
Dry friction, on the other hand, is ubiquitous in nature 
and plays a central role for diverse phenomena not only in 
physics and engineering, but also in biology and geology 
O In our model, the interplay of friction and noise 
leads not only to a non-zero mean velocity of the object 
in the absence of a mean force input, but also to non- 
monotonic and singular features in the transport proper- 
ties. By varying the dry friction strength three distinct 
modes of motion of the object can be induced: diffusive, 
directed and unidirectional. Strikingly, in the directed 
motion regime, the object's mean velocity increases when 
the dry friction coefficient increases, indicating that dis- 
sipation can in fact enhance transport. This highly coun- 
terintuitive effect might be exploited, e.g., in biological 
soft matter systems such as tethered vesicles or colloidal 
beads moving on lipid bilayers [TUlI12| . 

II. A MODEL FOR DIRECTED MOTION DUE 
TO DRY FRICTION AND SHOT NOISE 

The velocity v(t) of the solid object of mass m is 
described in the reference frame of the surface by the 
Langevin equation 

mv(t) = -jv - a(v)Am + £(t), (1) 

where £(i) denotes the random force. For simplicity, the 
discussion is reduced to one dimension. The friction be- 
tween the two solids is expressed phenomcnologically by 
two friction terms: a dynamic friction linear in v with 
strength 7, and a dry friction — a(v)Am with strength 
A. The sign function a(v), defined as a(v) = +1,0,-1 
for v > 0, ~ 0, < 0, respectively, is singular at v = 
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and models the non-analytic behaviour observed in the 
dynamics of a solid object due to dry friction. In par- 
ticular, the dry friction force in Eq. ([I]) can lead to com- 
plicated nonlinear stick-slip dynamics, where the object 
alternates in an oscillatory or chaotic way between stick- 
ing and sliding states [H |9] . We will sometimes find it 
useful to regularize the sign function by replacing it with 
a scaled hyperbolic tangent, where in the limit 

a(v) = lim tanh(u/e), (2) 

the singularity is recovered. Note that the strength of the 
dry friction is assumed directly proportional to the mass 
of the object, in accordance with the classical laws of 
friction dating back to the work by da Vinci, Amontons 
and Coulomb [S]. 

Models in the form of Eq. ([lj have been investigated re- 
cently both theoretically and experimentally for the cases 
when the vibrations £(i) are specified as asymmetric os- 
cillations [T3lU5| . or as Gaussian white noise with zero 
mean [T6Tf2"4] . In the first case, the dry friction effectively 
rectifies the asymmetric oscillatory vibrations, leading to 
a directed motion of the object. This rectification is due 
to the induced sticking of the object to the surface, so 
that a threshold force is needed to move the object. Di- 
rected motion results when the vibrations do not over- 
come the threshold symmetrically, even when the total 
force is zero on average. In the case of Gaussian noise, no 
directed motion results unless an additional bias, such as 
a constant force, acts on the object. As a generalization 
of Gaussian noise and in order to investigate the influ- 
ence of random spatial asymmetry, we consider Eq. ([I]) 
with taken as Poissonian shot noise (PSN). 

PSN is a mechanical random force, which is usually 
represented by a sequence of delta shaped pulses with 
random amplitudes A [7] 

C(t) = 2^(*-*fc)- (3) 

k=l 

The waiting time between successive pulses is assumed 
to be exponentially distributed with parameter A, where 
A is the rate at which pulses arrive. Then n*, the number 
of pulses in time t, follows a Poisson distribution 

P(n t =n) = ^£.e-* (4) 

with mean Xt. When a pulse occurs, its amplitude A is 
sampled, independently for each pulse, from a distribu- 
tion p(A). The mean and the covariance of £(i) are then 
given by 

«(*)) = MA), (5) 

(atw))-m))(m) = A(,4 2 >/(t-t'), (6) 

where (...)„ indicates an average over the amplitude dis- 
tribution p(A). In the following we assume that the noise 



does not exert a net force on the object, which requires 
a zero mean noise (£(t)) = 0. This can be achieved, e.g., 
by choosing a amplitude distribution with a vanishing 
mean. For an arbitrary p(A), the mean of the noise can 
be subtracted by hand, i.e., £(t) -> - (£(*))■ This 
case will be considered in more detail below. 

The PSN Eq. j3j converges to Gaussian white noise 
when A —> oo and (A) p — > with A^/l 2 )^ = const. 

[7]. In this limit, the model Eq. Qthus converges to 
Brownian motion with dry friction [16]. In the follow- 
ing, we focus on the non-Gaussian parameter regime of 
PSN and consider the stationary properties of the veloc- 
ity, which is sufficient for an understanding of the trans- 
port properties of the model. Even in the Gaussian case, 
the time-dependent statistics can only be determined ap- 
proximately using formal analytical methods 18-2T]. In 
Ref. [25] a similar dry friction model has been studied, 
where the noise is given by an asymmetric Markov pro- 
cess with zero mean that switches between two states. 
This type of noise also induces directed motion but no 
singular features arc reported. 

The velocity process Eq. (JlJ can be expressed more 
simply after division by m as 

v = -F(v)+£(t), (7) 

where the frictional forces are captured in 

F(v) = -v + a(v)A. (8) 

T 

The time scale r = m/7 is the inertial relaxation time. 
The noise £(t) now represents stochastic kicks with the di- 
mension of acceleration, which can formally be taken into 
account by an appropriate rescaling of the noise strength 
Aq — s- Ao/m. Taking the noise average in Eq. Q and 
considering the steady state where (v) = immediately 
leads to an expression for the mean velocity of the object: 

(v) = -Ar{a(v)) 

= At (^J p(v)dv - J p(v)dv\ . (9) 

Here, p(v) denotes the stationary distribution of the ve- 
locity, which has to be determined from the Kolmogorov- 
Feller equation associated with Eq. Q as discussed be- 
low. Eq. ([9| predicts that a positive mean velocity is in- 
duced when the total probability of observing a negative 
velocity is larger than the total probability of observing a 
positive one, which is somewhat counterintuitive, but is 
here a consequence of the particular form oip(v). More- 
over, Eq. Q is valid for generic noise sources with 
zero mean [31 and shows explicitly that a non-zero mean 
velocity of the object results from: 

(a) Inertia (non-zero r). 

(b) The influence of dry friction (non-zero A). 

(c) An asymmetric stationary velocity distribution. 
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Clearly, since the friction is symmetric, the asymmetry 
in the stationary distribution can only be induced by an 
asymmetric noise The asymmetric fluctuations then 
lead to directed motion of the object due to the presence 
of the nonlinear dry friction. For comparison, in an over- 
damped system, asymmetric PSN has also been shown 
to induce a macroscopic current in symmetric periodic 
potentials [25] . 

The velocity distribution p(v, t) of the velocity process 
Eq. Q can be derived from the Kolmogorov-Feller equa- 
tion associated with the PSN [13 [28] 



d_ 
at 



p(v,t) 



d_ 

dv 



F(v)p(v,t) 



+ X((p(v-A,t))-p(v,t)). (10) 



Under the conditions of stationarity and a zero probabil- 
ity current one obtains 



F(v)p(v) 



G(v - v')p(v')dv' 



(11) 



where the Green's function G(v) has the Fourier repre- 
sentation (indicated by a ~) 



G(k) = -iX 



( P{k) - 1 



(12) 



The integral representation Eq. (11) indicates that the 



PSN induces a non-local diffusion of the object. 



III. EXACT SOLUTION FOR ONE-SIDED 
SHOTS 

An analytical solution for the stationary distribution 
of Eq. ( 10 1 can be found when the amplitude distribution 



p(A) is given by an exponential distribution 



P(A) 



-A/A 



(13) 



where all amplitudes A are assumed to be positive, i.e., 
the PSN that we consider is one-sided. In order to obtain 
a noise with zero mean, we take for £ (t) 



Z(t) = Y^A k 5(t-t k )-\A , 



(14) 



fc=i 



so that (£(£)) = and the covariance is 2A/lQ<5(t — if). 
The noise thus consists of a random white shot noise part 
and a deterministic part that is equivalent to a constant 
negative drift force on the particle. Since the shot noise 
acts only in a one-sided fashion here, the noise £(£) is 
strongly asymmetric, even though its mean value is zero 
by construction. 

In the case of general values A and A the station- 
ary velocity distribution becomes strongly asymmetric, 



which is due to a lower cut-off for the possible velocity 
values at 



(A - A*)r 







< A < A* 



A > A* 



(15) 



where the critical dry friction strength A* is just the 
average of the PSN 



A* = A,4n 



(16) 



The lower bound in Eq. (151 can be understood as fol- 



lows. The stochastic kicks z(t) appear with rate A and 
only positive amplitudes. Between the kicks, the object 
will relax deterministically according to 



-v - a(v)A - A*, 
r 



(17) 



which represents the motion of an object under the in- 
fluence of dry friction and a constant force —A*. This 
equation has two different fixed points depending on A. 
For A < A*, the fixed point is at v* = (A — A*)t, which 
follows from setting v — in Eq. (17) and solving for v. 
However, for A > A* there is always a net force oppo- 
site to the direction of motion and the fixed point is at 
v* = 0. Physically, the object will not be able to move 
with a negative velocity in this regime because the dry 
friction dominates the constant force and the object will 
always get stuck over time when relaxing from a stochas- 
tic kick. 

Both regimes result in a non-zero mean velocity as will 
be discussed in the following. To summarize at this point, 
there are three distinct modes of motion of the object: 

1. Diffusive motion ((v) — 0) when A = 0. 

2. Directed motion ((v) > 0) when < A < A*. 

3. Unidirectional motion ((«) > and v(t) > 0) when 
A > A*. 

Diffusive motion also occurs for non-zero A in the Gaus- 
sian limit of the PSN, which has been discussed in 
Refs. [I"8ll21j . For A and Aq held constant, so that 
A* = const., the transitions between the three regimes, 
and hence from diffusive to unidirectional motion, take 
place upon increasing the friction coefficient A. In other 
words, the transport on the surface in fact improves for 
increasing dry friction. In all three regimes of motion, in- 
dividual velocity trajectories reveal transitions between 
sticking and sliding states for non-zero A, which is a char- 
acteristic feature of motion subject to dry friction (figures 
not shown). In the unidirectional motion regime, Eq. @ 
for the mean velocity must break down, since it would 
predict a negative mean even though the object's instan- 
taneous velocity is never negative. The reason for the 
breakdown is the appearance of a delta-peak in p(v) at 
v = 0. 
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To explore these effects quantitatively, we consider the 



Kolmogorov- Feller equation (10), which under our as- 



sumption of one-sided exponential shot noise becomes 
8 



9 , , 



dv 
-A 



[F(v)+A*]p(v,t) 
d 1 



dv 1 + A d/dv 



p[v,t). (18) 



Interestingly, the stationary state of this equation can 
be found exactly for arbitrary F(v) using an operator 
inversion, leading to p(v) oc h(v) with |29j 



h(v) 



X 



F(v) 



exp 



v 



X 



F(v') + A* 



dv' }(19) 



The stationary distribution can therefore be determined 
not only for the piecewise linear friction Eq. ^ , but also 
for the smooth representation Eq. ^ of the a(v) singu- 
larity. This allows for a validation of the singular effects 
of the dry friction by a limit procedure. For F(v) given by 
Eq. Q, the exact expression for the integral in Eq. (19) 
can be determined separately for v > and v < 0. The 
solution is then 



h(v) = 



{ x 






A*+A 




1 » 






I. A*-A 





Ar-1 



e A " , v > 



At- 1 



(20) 



e A o , u_ < w < 0, 



where u_ is given in Eq. ( 15 ) and v. 



-(A + A*)t (so 



that v+ is the hypothetical fixed point of Eq. (17) that 
one would obtain by replacing a(v) by 1). For A = 0, 
i.e., in regime 1, p(v) is identical for v > and v < as 
expected, because there is then no singular effect from 
the dry friction. 



IV. SINGULAR FEATURES 



Eq. (20) is plotted in Fig. [T] together with results from 
a direct simulation of the equation of motion (JlJ. The 
simulation uses a PSN increment method developed in 
and the smooth tanh representation Eq. |2]). For 
comparison, smooth versions of p(y) are also shown in 
the inset. One clearly recognizes the trend towards a 
discontinuity at v — as e — » 0. The gap at v = 
satisfies 



p(Q- 



A 



p(0H 



A' 



(21) 



and increases monotonically with A in the interval < 
A < A*. As A — > A* the gap becomes infinite while 
V- — » 0, i.e., the support of p(v) for v < vanishes. This 
indicates that the distribution exhibits a delta peak at 
v = for A = A*. This peak in fact persists for A > A* 
as numerics show and can be obtained also from Eq. ( 19 1, 
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FIG. 1: (Colors online) The normalized stationary velocity 
distribution Eq. ( 20 1 for two different values of A in the regime 
< A < A* . One notices the discontinuity in the distribution 
at v = 0. The dots represent simulation results. Inset: the 
stationary velocity distribution using the tanh representation 
of the a(v) term, Eq. pL and A = 4. Parameter values: 
A = 5, A Q = 5, r = 1. 



For A > A* one can show that the lower cut-off V- on v 
in such a regularized model is negative and of 0(e). The 
term — v/Aq in the exponent in Eq. (19) is negligible be- 
tween this cut-off and v = 0, and the total mass of h(v) 
in this rang e becomes l-exp{- / _ X/(F(v') + A*)dv'}. 
The second term here vanishes because the integral is 
logarithmically divergent at the lower end. The distribu- 
tion p(v) thus contains finite probability mass within an 
O(e) range of v = 0, and this becomes a delta-peak for 
e 0. 

The stationary solutions of diffusion processes are usu- 
ally assumed to be continuous, which is required by the 
local nature of ordinary diffusion. In our case, PSN in- 
duces a non-local diffusion of the object, so that continu- 
ity of the solution is not required. For comparison, in the 
limit of white Gaussian noise £ (t) the Kolmogorov-Feller 
equation (10) reduces to the standard Fokker-Planck 



equation and Eq 
to p'(v) on the rig 



11 ) would contain a term proportional 
In this limit the station- 



;ht-hand side. 

ary distribution can be found in a straightforward way 



p(v) oc exp 













-Ha)} 









(22) 



for a noise strength D = XA$. Eq. (22) is continuous, 
but exhibits a cusp singularity at v = 0. Due to the 
symmetry of p(v) the mean velocity is zero in this case 
and the transport is purely diffusive. 

In order to understand quantitatively the gap in p(v) 
as well as the appearance of the delta peak for A > A* we 



by considering a regularized version of cr(y) as follows. consider the integral equation ( 11 ). For the noise Eq. (14) 
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with exponentially distributed amplitudes, the Green's 
function is given by G(v) — Q(v)\e~ v / A ° . Eq. (11) thus 
reads 

(F(v) + A*)p(v) = X [ e-( v - v 'V Ao p(v')dv'. (23) 



Assuming first that p(v) consists of two separate parts 
for v > and v < 0, but no delta-peak at v = 0, the 
right-hand side is continuous and so the values of p(v) 
either side of the gap at v = satisfy the relation 



(A + A*)p(0+) = (A* - A)p(0" 



(24) 



The ratio of the gap Eq. (21 1 follows immediately, if < 
A < A*. For A > A* the cut-off at v*_ — enforces 
p(0~) = so that Eq. d24 1 would require p(0 + ) = 0. 
However, a zero value of p{v) at v — + is not consistent 
with the solution for v > in Eq. (pOl. The conclusion is 



that there has to be an additional contribution in order 
to satisfy Eq. (24), which can only come from a delta 
peak in p{v) at v = 0. Therefore, we assume that, for 
A > A* (regime 3), the stationary distribution has the 
form 



p(v)=T S(v) + e(v)c + h(v), 



(25) 



where c+ is a normalization constant. Eq. (24) is then 
replaced by 



(A + A*)c+/i(0 H 



Ar , 



(26) 



so that together with the normalization condition 
J p{v)dv — 1, the two unknowns r and c + can be 
uniquely determined. This yields c + = Tq and Tq — 
(1 + J °° /i(w)dw) -1 . The distribution obtained in this 
way agrees well with simulation results (cf. inset of 
Fig. [2]) . The transition from the distribution Eq. ( 20 ) 
for < A < A* to Eq. (l25|) for A > A* occurs continu- 



ously, since the delta-peak amplitude Tq is just the area 
of the v < part of p(v) in Eq. (jp) as A ->■ A*. This 

follows from the limit: limA->.A* /„ h(v)dv — 1. 

In the case A > A*, Eq. Q is not valid, but an anal- 
ogous simple expression for the mean velocity can be de- 
rived using a limit procedure. Using the smooth tanh 
representation Eq. |2]), the mean velocity of the process 
Eq. |7| can be written in the form 



(v) 



At / tanh 



Q h(v)di 



At I tanh(-) h(v)dv, (27) 



where w_ is the lower cut-off as before. The second in- 
tegral becomes c + h(v)dv as e — > + . For the first 
integral, one can show that 



lim 



tanh 



(^) h{v)dv 



= -A*Tn 



(28) 




FIG. 2: (Colors online) Plot of the mean velocity (v), given by 
Eqs. (9]) and |29| for < A < A* and A > A*, respectively, 
as a function of A/A* for three different values of a = Ar. At 
A = A* the crossover from the directed to the unidirectional 
regimes of motion leads to a cusp singularity. Inset: Plot 
of the stationary velocity distribution Eq. ( |25[ ) for A > A*. 
The peak at v — stems from the delta-peak with amplitude 
To. The dots represent simulation results. Parameter values: 
A = 1.5A*, A = 4, A = 1, t = 1. 



when A > A*. The analogue of Eq. Q is then 

/>oo 

(v) = A*rr - Arc+ / h(v)dv. (29) 
Jo+ 

In order to investigate the properties of the mean veloc- 
ity we rewrite the expressions for (v), Eqs. ^ and (29), 
in terms of the three parameters Aq, A/A*, and a = At. 
This leads to the form (v) — AoG(-^,a), where G has 
a different functional form for < A < A* and A > A*. 
The mean velocity rescaled by Aq as a function of A /A* 
is plotted in Fig. [2] for three different a values and shows 
a non-monotonic behavior. The crossover between the 
directed and unidirectional regimes of motion leads to a 
cusp singularity in (v) at A = A*. In the directed mo- 
tion regime, (v) increases with A up to a maximal value 
of (v) — Aq as a — > 00, indicating that the transport gen- 
erally improves on a rougher surface for the same noise 
characteristics (fixed A and Aq) in this regime. Opti- 
mal transport is achieved for A = A* and a — > 00. In 
other words, the object attains the maximal mean veloc- 
ity when the mean of the noise balances the dry friction 
force (for v > 0) and the object has no time to relax in 
between successive noise pulses. Eventually, in the uni- 
directional motion regime, (v) decays to zero as A — > 00. 



V. CONCLUSION 

We have investigated an exactly solvable nonlinear 
model for the directed motion of an object due to dry 



G 



friction and shot-noise with zero mean. The transport in 
this model behaves in a non-monotonic way as the dry 
friction strength is increased, exhibiting a transition to 
a unidirectional mode of motion at a critical dry friction 
strength A*. This transition is indicated by a cusp sin- 
gularity in the mean velocity. The appearance of such 
a singularity at the transition between two qualitatively 
different modes of motion is reminiscent of a dynami- 
cal phase transition, which occurs here in a simple one- 
dimensional transport model. 

A large variety of different distributions for the ampli- 
tudes of the shot noise can be incorporated in our model 
Eq. 0. In the presence of negative amplitudes, which 
can be realized, e.g., by choosing a shifted Gaussian or 
exponential amplitude distribution, the stationary dis- 
tribution of the velocity could not be obtained in closed 
analytical form so far. However, some of the distinct sin- 
gular features discussed here, namely the discontinuity 
and the delta peak in p(v) at v = are a general con- 
sequence of the dry friction nonlinearity coupled to non- 
local shot noise. For this type of noise a path-integral 
solution similar to the one for Gaussian noise, derived 
in Refs. [TH E], can also be developed. This will be 
discussed in a forthcoming article. 

The results presented could be tested experimentally 
in setups similar to the ones used by Chaudhury et al 



[22TI24) . Instead of Gaussian white noise vibrations one 
would have to induce vibrations with PSN statistics. By 
balancing the one-sided noise with a constant force, e.g., 
gravitation, a zero mean of the resulting force could be 
achieved. It would be interesting to observe in particular 
the sharp transition between the directed and unidirec- 
tional motion regimes. Moreover, the enhancement of 
transport for increased dissipation might be exploited in 
artificial or biological soft matter systems. Recent exper- 
imental results on diffusion statistics of tethered vesicles 
[TOl ITT] or colloidal beads moving on lipid bilayers [T2] 
have been shown to agree with an effective dry friction 
model in the form of Eq. ^ [21 . The PSN used in our 
model might be realized in these systems by an addi- 
tional active chemical process that occurs on a separate 
time scale. To what extent the properties of our model 
persist in the presence of thermal noise remains to be 
investigated, but since thermal noise does not impose a 
force bias the qualitative features of the mean velocity 
should agree with our results. 
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